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Each  year  the  Aeronautical  Systems  Division  (ASD),  Air  Force  Systems  Command 
(AFSC),  sponsors  a  Science  and  Engineering  Symposium  in  advance  of  the  Annual  Air 
Force  Science  and  Engineering  Symposium.  This  provides  a  specific  motivation  for  ASD 
personnel  to  prepare  papers  that  reflect  the  results  of  their  efforts.  The  variety  of  sub¬ 
jects  also  provides  an  opportunity  for  interdisciplinary  exchange  of  information  that  is 
becoming  ever  more  important. 

This  year  the  symposium  papers  are  being  published  individually  to  facilitate  distribu¬ 
tion  and  retention.  However,  each  paper  carries  this  same  foreword  which  lists  the  titles 
of  all  papers  together  with  the  authors  and  the  ASD  Technical  Documentary  Report  (TDR) 
numbers.  Readers  who  are  interested  in  obtaining  copies  of  other  papers  are  urged  to 
contact  the  authors  directly  or  the  Defense  Documentation  Center,  Alexandria,  Virginia. 

It  should  be  noted  that  certain  papers  are  classified  and  are  available  to  only  those  per¬ 
sons  having  proper  security  clearances  and  a  “need-to-know.” 

This  paper  is  one  of  21  presented  at  the  "ASD  1963  Science  and  Engineering  Symposium” 
held  at  Wright-Patterson  Air  Force  Base,  Ohio,  18-19  September  1963.  They  consist  of  17 
CONTRIBUTED  and  4  INVITED  papers,  listed  below.  *The  5  contributed  papers  that  are 
asterisked  were  also  presented  at  the  10th  Annual  Air  Force  Science  and  Engineering 
Symposium  held  at  the  Air  Force  Academy,  Colorado  Springs,  Colorado  on  8,  9  and  10 
October  1963. 
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ABSTRACT 


This  study  considers  the  concept  of  a  hypersonic  glider-type  spacecraft  utilizing  its 
aerodynamic  maneuvering  capability  in  performing  orbital  plane  change.  For  lifting 
vehicles  an  optimization  procedure  is  developed  which  defines  the  proper  vehicle  attitude, 
propulsion  utilization  and  sequence  of  operations  to  produce  the  maximum  plane  change 
for  a  given  fuel  expenditure.  The  results  obtained  are  compared  with  the  fuel  requirements 
for  a  pure  propulsion  (nonlifting)  plane  change  while  remaining  in  orbit.  Specifically,  the 
optimum  bank  angle,  angle  of  attack,  entry  angle,  thrust  alignment  and  thrusting  procedures 
are  defined.  In  addition,  the  advantages  of  high  L/D  vehicles  are  graphically  illustrated. 

The  method  is  seen  to  be  more  efficient  than  the  pure  propulsion  method,  but  is  found  to  be 
far  more  complex  and  requires  longer  times  to  execute. 
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I.  INTRODUCTION 


During  space  flights  of  extended  duration,  or  to  fulfill  specified  mission  objectives,  it 
will  become  necessary  for  future  spacecraft  to  possess  an  orbital  plane  change  capability. 

In  general,  there  are  two  ways  in  which  this  may  be  accomplished:  (1)  by  a  pure  trans¬ 
verse  propulsive  thrust  while  remaining  in  orbit,  or  (2)  by  atmospheric  maneuvering  uti¬ 
lizing  aerodynamic  lift  and  small  amounts  of  selected  propulsion  after  descent  from  orbit, 
followed  by  reascent  to  orbital  altitude.  The  first  method  is  relatively  simple  in  execution, 
but  expensive  in  terms  of  fuel  consumption.  The  second  method,  titled  “synergetic  or  aero- 
propulsive  maneuvering,”  can  be  more  economical  in  its  fuel  requirements,  but  is  far  more 
complex  from  a  guidance  and  control  standpoint  and  requires  longer  durations  to  execute. 

While  utilization  of  the  pure  propulsion  method  results  in  nearly  instantaneous  plane 
change,  the  amount  of  propellant  needed  may  be  restrictively  high.  For  example,  assuming 
a  circular  orbit  of  100  nautical  miles  altitude,  a  transverse  velocity  impulse  (which  changes 
the  direction  but  not  the  magnitude  of  the  original  velocity  vector)  of  approximately  447  ft/ 
sec  is  required  for  each  degree  change  in  the  orbital  plane  angle.  While  this  may  not  appear 
excessive  at  first  glance,  consider  that  this  means  velocity  impulses  of  4460  ft/sec  for  a 
10-degree  plane  change  and  13,243  ft/ sec  for  a  30-degree  plane  change.  In  terms  of  pro¬ 
pellant  consumption  for  a  maneuvering  system  having  an  1  =  310  seconds,  0.564  pounds 

sp 

of  fuel  per  pound  of  inert  structure  is  required  for  a  10-degree  plane  change  and  2.77 
pounds  of  fuel  per  pound  of  inert  structure  is  required  for  a  30-degree  plane  change.  While 
plane  changes  of  the  magnitudes  cited  above  may  not  be  required  in  all  space  missions,  the 
above  consideration  of  the  magnitudes  of  impulse  and  propellant  consumption  yields  an  idea 
of  the  problem  involved. 

The  method  utilizing  atmospheric  maneuvering  in  conjunction  with  propulsion  augmenta¬ 
tion,  termed  the  “synergetic”  method  (Ref.  1),  is  the  subject  of  this  report.  Due  to  the 
complexity  of  the  maneuver,  it  must  be  dealt  with  in  several  distinct  phases,  and  is  only 
considered  because,  properly  employed,  it  is  more  economical  in  terms  of  propellant 
consumption,  a  factor  of  primary  importance  in  present-day  technology.  This  advantage 
is  derived  from  the  fact  that  for  near-earth  orbits  it  is  cheaper  to  alter  altitude  than 
direction. 

The  analysis  is  performed  in  separate  and  distinct  steps,  and  throughout  an  effort  is 
made  to  optimize  the  procedure  so  that  the  minimum- energy  path  is  followed  as  closely 
as  possible.  The  individual  steps  comprising  the  method  are  (Fig.  1): 

(1)  retrograde  impulse  and  deorbit; 

(2)  reentry; 

(3)  circularization  of  the  orbit  while  forcing  the  velocity  vector  to  y  =  0°; 

(4)  thrusting,  banked  turn  within  the  atmosphere  during  which  the  major  portion  of  the 
plane  change  is  accomplished; 

(5)  velocity  impulse  sufficient  to  exit  the  atmosphere  and  return  to  the  original  orbital 
altitude;  and 
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(6)  velocity  impulse  to  circularize  and  reestablish  the  orbit. 

It  is  the  purpose  of  this  report  to  ascertain  whether  the  added  complexity  and  time  re¬ 
quirement  of  this  method  can  be  justified  in  terms  of  the  quantity  of  propellant  saved. 


Figure  1.  Schematic  of  Synergetic  Maneuver 
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11.  ANALYSIS 


Before  starting  the  actual  analysis,  it  is  necessary  to  list  certain  assumptions  which 
govern  the  entire  investigation.  It  is  assumed  first,  that  an  accurately  throttleable,  multi¬ 
start-capability  rocket  engine  is  available;  second,  that  thrust  applications  in  space,  to 
circularize  the  orbit  after  reentry  and  to  exit  the  atmosphere,  result  in  instantaneous 
velocity  and/or  direction  changes;  third,  that  the  vehicle  is  stable  (ac  =  eg)  so  that  no 
allowance  is  necessary  for  thrust  vector  control;  and  fourth,  the  chosen  L/D  ratios  are 
attainable  by  present  or  proposed  spacecraft.  Since  each  step  in  the  maneuver  is  a  sep¬ 
arate  entity,  each  is  treated  as  such  and  the  resultant  velocity  increments  summed  to 
yield  the  total  propulsion  requirements  for  a  given  plane  change  magnitude.  In  an  effort  to 
present  as  thorough  a  study  as  possible,  covering  the  greatest  variety  of  circumstances, 
the  following  conditions  are  assumed  throughout  the  report: 

(1)  the  original  orbital  altitude  is  100  nautical  miles  since  Reference  1  shows  that  lower 
orbital  altitudes  are  more  efficient  for  synergetic  maneuvering  and  this  is  a  reason¬ 
able  minimum  from  an  orbital  lifetime  standpoint; 

(2)  the  aerodynamic  characteristics  may  be  expressed  by  Newtonian  techniques  (Appen¬ 
dix  I); 

(3)  the  original  vehicle  weight  and  reference  area  remain  the  same  for  all  L/D  values 
since  variations  in  these  parameters  result  in  small  performance  changes  as  shown 
in  Appendix  II; 

(4)  no  thermal  or  structural  limitations  exist.  (However,  maximum  temperatures  during 
the  maneuver  are  presented  for  comparative  purposes.) 


DLORB1T  PHASE 

From  a  100- nautical  mile  circular  orbit,  a  deorbit  is  required  such  that  the  vehicle 
reaches  300,000  feet  at  a  specified  reentry  angle  with  minimum  energy  loss.  The  minimum- 
energy-loss  transfer  between  coplanar  orbits  has  been  shown  (Ref.  2)  to  be  a  Hohmann 
transfer  ellipse  with  apogee  and  perigee  tangent  to  the  original  and  final  orbits  respectively. 
The  retrograde  velocity  necessary,  and  the  resultant  entry  velocity  for  this  maneuver  can 
be  calculated  from  the  Keplerian  equations  (Ref.  3)  with  negligible  error.  From  the  con¬ 
stancy  of  energy  and  angular  momentum  at  any  point  on  a  given  orbit,  one  obtains 


and 


r2  v2 


=  r  V  cos  y 
en  en  h 


en 


(1) 

(2) 


Combining  these  equations,  recognizing  that  rs  =  ra ,  and  solving  for  the  entry  velocity, 

V  ,  obtain 
en 
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where 


Using  Equation  (3),  the  entry  velocity  at  300,000  feet,  after  descent  from  a  100-nautical 
mile  circular  orbit,  may  be  calculated  for  the  specified  entry  conditions.  Then,  from 
Equation  (2)  the  value  of  V2  is  obtained.  Thus,  the  retrograde  impulse  required  is,  by 
definition 


AV  =  V,  -  V2  .  (4) 

To  ascertain  that  this  is  the  minimum  A V  value,  consider  the  following  sketch: 


From  this  sketch 


V2  cos  y  =  V,  +  Av  cosX.  . 
For  any  orbit,  the  constancy  of  angular  momentum  yields 

r  V  cosy  =  Constant  . 


(5) 

(6) 


Therefore 

rz  v2  C0SX  =  ren  ven  cosren 
Combining  Equations  (5)  and  (7)  obtain 

7T  Ven  Ten 
cos  X 

Since  the  conditions  at  the  time  of  retrofire  are  fixed  by  the  original  orbit  and  those  at 
atmospheric  entry  are  established  by  definition,  Equation  (8)  may  be  rewritten  as 


A  Constant 

Av  =  — — - 

cosX 


(9) 
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Thus,  if  IavI  is  to  be  a  minimum  for  a  given  reentry,  then 


cosX  must  be  a  maximum,  or 


a  pure  retrograde  impulse. 


X  =  180°  , 


(10) 


The  value  of  the  retrograde  impulse  obtained  by  Equations  (3)  and  (4)  is  shown  in  Figure 
2.  These  values,  then,  are  the  minimum  velocity  impulses  necessary  to  perform  the  de- 


0  0.5  1.0  1.5 

-Xen  (deg) 

Figure  2.  Variation  of  Retrovelocity  with  Entry  Angle  (100-NM  Orbit) 


REENTRY  PHASE 

At  an  altitude  of  300,000  feet  aerodynamic  effects  become  appreciable  and  it  is  neces¬ 
sary  to  use  more  refined  techniques  in  the  analysis.  It  is  also  desirable  to  know  what 
reentry  procedures  will  lead  to  minimum  energy  dissipation  for  a  given  plane  change;  that 
is,  what  bank  angle,  entry  angle,  angle  of  attack  and  pullout  procedure  will  require  the 
least  amount  of  propulsion  to  attain  the  same  results.  To  take  all  of  the  possibilities  into 
account,  computer  programs  were  set  up  based  on  the  entry  angles  and  velocities  obtained 
from  the  DEORBIT  PHASE  analysis.  Cases  were  run  for  several  reentry  angles  from  -0.1° 
to -1.5°,  at  various  bank  angles  from  0°  to  90°  and  for  various  angles  of  attack  based  on 
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(L/D)  values  from  1  to  4.  The  aerodynamic  characteristics  for  the  chosen  (L/D) 

'  '  'max  1  '  '  'max 

values  were  computed  by  means  of  Newtonian  aerodynamics  as  presented  in  Appendix  1. 

The  results  of  the  computer  runs  defined  the  optimum  angle  of  attack  and  bank  angle  for 

the  reentry  phase  as  shown  in  Figures  3  and  4  respectively.  A  comparison  of  Figures  3 

and  1-1  shows  that  the  optimum  angle  of  attack  for  the  reentry  phase  is  that  for  (L/D)  , 


0  I  2  L  3  A 

~D~ 

Figure  3.  Optimum  Angle  of  Attack  Versus  L/D 


L_ 
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Figure  4.  Optimum  Bank  Angle  Versus  L/D 
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The  above  reentry  conditions  were  run  to  pullout,  or  y  =  0°.  After  considering  this  data 
it  was  felt  a  more  economical  pullout  maneuver  could  be  effected  in  light  of  the  magnitude 
of  the  velocity  impulse  required  before  the  next  phase  of  the  maneuver  could  be  accomp¬ 
lished.  At  selected  points  during  reentry,  the  velocity  impulse  was  computed  which  placed 
the  vehicle  at  circular  speed  and  y  =  0°  by  using  the  law  of  cosines, 

Av  =  /v2  +  Vc2  -2 V  vc  cosy  (11) 

where 


V  =  local  inertial  velocity 
Vc  =  local  circular  velocity 
y  =  local  flight  path  angle. 

These  velocity  impulses  were  calculated,  and  the  minimum  determined,  thus  defining  the 

pullout  altitude  and  the  starting  point  for  the  next  phase  of  the  analysis.  Figure  5  presents 

the  pullout  flight  path  angle  and  Figure  6  is  a  plot  of  pullout  altitude  versus  entry  angle 

for  the  various  L/D  values  for  0  and  a  .  These  values  were  defined  by  the  point  of 
'  opt  opt  3  1 

minimum  velocity  impulse  for  circularization  which  is  shown  in  Figure  7.  A  point  of  note 
in  Figure  7  is  the  plotting  of  negative  velocity  impulses.  These  are  recognized  as  improper, 
and  are  merely  a  result  of  interpretation.  The  cases  in  which  they  occur  reached  y  =  0° 
aerodynamically  and  hence  the  impulse  served  only  to  slow  the  vehicle  to  circular  speed. 

In  practice  this  is  not  done,  but  the  turn  is  performed  at  supercircular  speed  and  the  excess 
serves  to  reduce  the  impulse  needed  to  exit  the  atmosphere.  The  heading  change  accomp¬ 
lished  during  reentry  is  shown  in  Figure  8.  After  pullout  is  completed,  the  major  portion 
of  the  plane  change  must  be  undertaken. 


ATMOSPHERIC  MANEUVERING  PHASE 

After  pullout,  the  atmospheric  maneuvering  begins  in  order  to  accomplish  the  major  por¬ 
tion  of  the  plane  change.  As  suggested  in  Reference  1,  thrust  is  maintained  equal  to  drag 
during  the  atmospheric  maneuvering  so  that  no  velocity  loss  is  incurred.  This  procedure 
results  in  a  thrusting  banked  turn  at  constant  velocity,  angle  of  attack,  flight  path  angle  and 
bank  angle. 
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Considering  the  above  sketch,  the  basic  equations  governing  the  aerodynamic  maneuver 


V  \f/  cos  y  =  L  sin  <jfc>  +  T  sinA  sin<£ 


i  Vy 


mV  cosy 


+  Lcos(£-mg  cosy  +  T  sinA  coscp 


(i: 


mV  =  -D  —mg  siny  +  T  cosA 


m  Vex  -m  Veffectjve 


(1' 

(1! 


Figure  5.  y  for  Minimum  AV  Versus  Entry  y  for  Various  L/D  Values 
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Employing  the  identity 


d  [f/ 


d  f  dm 


dt  dm  dt 

and  using  Equation  (15),  Equation  (17)  can  be  written,  after  some  manipulation  as 


d  \j/ 


—  d  ( J^n  m) 


=  V 


ex 


cosA  +  s  i  n A  j 


iin<£ 


(20) 


(21) 


Tlie  greatest  change  in  i(i  for  minimum  fuel  expenditure  may  be  obtained  by  differenti¬ 
ating  Equation  (21)  with  respect  to  the  independent  variables.  Differentiation  with  respect 
to  A  produces  the  following  optimum: 


cot  A„_ .  =  — 

opt  o 


(22) 
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Figure  7.  Minimum  AV  Required  to  Circularize  Orbit  During  Reentry  for  Various 
L/D  Values 
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Figure  8.  Change  in  Heading  During  Reentry  for  Various  L/D  Values 
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Using  Equation  (22),  Equation  (21)  becomes 


d  \j/ 

- d  iJn  m  ) 


si  ncj> 
V 


(23) 


It  is  now  apparent  that  the  most  efficient  turning  rate  may  be  achieved  by  flying  at 
(L/D)max  or  .  The  optimum  values  of  a  and  A  are  listed  in  Table  1  and  in  Figures 

1-1  and  9  respectively,  for  various  L/D  values. 


TABLE  1 


OPTIMUM  A  VALUES 


L 

D 

aopt 

(degrees) 

^opt 
(degrees ) 

■^opt  aopt 

(degrees) 

1 

29.25 

45.00 

15.75 

1.5 

22.17 

33.67 

1L.50 

2 

17.51 

26.57 

9.06 

2.5 

14.35 

21.8 

7.45 

3 

12.20 

18.47 

6.27 

3.5 

10.6 

15.95 

5.35 

4 

9.35 

14.03 

4.68 

Thus  for  maximum  performance,  the  engine  should  be  placed  in  the  vehicle  declined 

from  the  vehicle  centerline  at  an  angle  of  (A  -  a  ). 

b  v  opt  opr 

Referring  to  Equation  (23),  it  is  noted  that  sin  0/V  must  also  be  maximized  for  peak 
efficiency.  Equations  (18)  and  (19)  may  be  used  to  eliminate  V  from  Equation  (23).  Thus 
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L_ 

D 


Figure  9.  Optimum  Thrusting  Angle  Versus  L/D 


Differentiation  of  Equation  (24)  with  respect  to  0  produces 

I 


COS 


opt 


3B 


+  V  (lB  )  + 


(25) 


From  Equation  (25)  the  variation  in  0  with  B  was  calculated  and  the  results  presented 
in  Figure  10.  The  dependency  of  0  on  altitude  and  L/D  was  also  computed  and  plotted  in 
Figure  11.  Thus,  substitution  of  0  from  Figure  11  into  Equation  (24)  results  in  the  maxi¬ 
mum  turn  for  a  given  fuel  expenditure.  Utilizing  Equation  (16),  Equation  (24)  becomes 

sin<£ 


Ai// 


Av 


ef  f 


it)  +  ' 

k  ^  'opt 


opt 


['+  B  cos<£0pt  ] 


(26) 


which  is  plotted  in  Figure  12. 


EXIT  PHASE 

After  a  prescribed  heading  change  has  been  achieved,  a  velocity  impulse  must  be  added 
sufficient  in  magnitude  to  propel  the  vehicle  from  the  atmosphere  on  an  elliptical  path  hav¬ 
ing  an  apogee  at  100  nautical  miles.  To  maintain  the  concept  of  minimum  energy  expended, 
the  vehicle  must  exit  the  atmosphere  on  a  Hohmann  transfer  ellipse.  The  optimum  bank 
angle  and  angle  of  attack  for  reentry  were  also  adopted  for  the  ascent  phase.  The  velocity 
impulse  required  to  insert  the  vehicle  on  the  proper  Hohmann  transfer  ellipse  after  complc' 
tion  of  atmospheric  maneuvering  could  not  be  determined  by  a  direct  method. 
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Figure  10,  Optimum  Bank  Angle  Defined  by  Flight  Conditions  for  Atmospheric 
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Figure  11.  Optimum  Bank  Angle  Versus  Pullout  Altitude  for  Various  L/D  Values 
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Figure  12.  Plane  Change  Efficiency  During  Thrusting  Banked  Turn  for  Various 
L/D  Values 

Therefore,  a  solution  was  obtained  by  an  iteration  technique  in  which  the  exit  trajectory 
and  Hohmann  transfer  ellipse  equations  were  solved  simultaneously.  The  results  of  this 
analysis  in  terms  of  velocity  requirements  are  plotted  in  Figure  13.  The  heading  change 
for  the  exit  phase  is  shown  in  Figure  14. 


ORBIT  REESTABLISHMENT  PHASE 

Once  the  proper  exit  conditions  are  obtained  it  is  a  simple  task  to  calculate  the  apogee 
velocity  of  the  resultant  ellipse  using  Equation  (2)  subscripted  for  exit  conditions.  The 
difference  between  this  calculated  value  and  the  known  circular  velocity,  then,  is  the  re¬ 
quired  recircularization  velocity  impulse  which  completes  the  maneuver.  This  velocity 
increment  is  given  in  Figure  15. 

COMBINATION  OF  PURE  PROPULSION  AND  AERODYNAMIC  1URN 

Throughout  the  preceding  analyses,  the  most  efficient  turning  path  was  sought  and  fol¬ 
lowed.  Consider  now  the  application  of  additional  thrust  normal  to  the  orbital  plane,  thus 
obtaining  small  increments  of  A 0  while  deorbiting,  circularizing,  ascending  and  reestab¬ 
lishing  orbit.  The  advantage  is  readily  apparent  if  one  considers  a  30°  right  triangle  in 
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which  a  transverse  AV  of  0.500  is  obtained  for  a  modest  reduction  of  the  longitudinal  AV 
to  0.866  for  the  same  total  velocity  increment. 

Keeping  the  above  in  mind,  an  optimization  of  this  principle  is  attempted.  Consider  the 
following  sketch. 


Figure  13.  Velocity  Impulse  Required  to  Exit  Atmosphere  and  Reach  100-NM  Apogee 
for  Various  L/D  Values 
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Figure  15.  Velocity  Impulse  Required  to  Reestablish  the  Orbit  (100  NM)  for  Various 
L/D  Values 


From  V) ,  a  new  velocity  Vn  is  to  be  obtained  such  that  Ai/j/AV  is  a  maximum.  Thus,  from 
the  preceding  sketch 


and 


V2  sin  Ax//  =  AV  sinS 
V,  +- AV  cosS  =  V2  cos  A  v// 


For  an  optimum  maneuver, 


(— ) 

'  AV  'total 


maximum  . 


(27) 


(28) 


(29) 


Thus,  (Ai j) /AV)  is  composed  of  the  impulse  terms  and  the  aerodynamic  turning 

.  total 

term  or 
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(—) 

\  Av/  h 


,  A vK  /  Avi/  v 

("Av" ) imp  AV'nnp  +  (~Av), 


AV, 


aero 


aero 


Av  'total 
Then,  from  Equations  (27)  and  (26) 


(30) 


AV.  +  AV 
imp  aero 


a  i  si n  S 

A\//  «  -  —  AV-  n  + 

T  \j9  imp 


y^f 


Av 


aero 


(31) 


where,  in  Equation  (26)  the  value  of  0^  is  taken  as  90°  and  in  Equation  (27)  small  angles 
of  A 0  are  assumed.  Also,  from  Equation  (28) 

v2-v, 


Av.  « 

imp 


cos< 


(32) 


and  by  definition 


Av 


aero 


=  AV 


total 


Vg  ~V| 
cos8 


Using  Equations  (31)  and  (32),  Equation  (30)  may  be  written  as 


(33) 


A  \f/  =  t  a  n  8 


V2  -Vi 
V, 


(34) 


Maximizing  A0  for  a  specific  AV  with  respect  to 


yields, 


where  V2  ^  V  has  been  assumed, 
c 

Thus,  using  Equation  (35),  the  optimum  angle  to  offset  the  thrust  is  defined  and  these 
values  are  shown  in  Figure  16.  It  was  found,  because  of  the  complex  orbital  parameter 
interactions,  that  this  offset  may  be  efficiently  used  only  during  pullout  and  exit  and  not 
during  deorbit  and  reestablishment  as  was  first  assumed.  The  velocity  requirement  and 
the  heading  change  obtained  by  thrust  offset  at  pullout  and  exit  are  shown  in  Figures  17 
and  18  respectively. 
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III.  DISCUSSION  OF  RESULTS 


At  this  point  the  conclusions  reached  in  the  Analysis  Section  on  the  most  efficient  tech¬ 
nique  for  accomplishing  the  orbital  plane  change  are  summarized.  As  pointed  out  in  that 
section,  a  pure  retrovelocity  impulse  yields  the  most  economical  reentry  in  terms  of  pro¬ 
pellant  consumption.  This  lends  further  support  to  the  conclusion  that  it  is  easier  to  alter 
altitude  than  direction  while  in  near-earth  orbits.  As  reentry  begins,  however,  it  is  neces¬ 
sary  to  consider  the  effects  of  entry  angle  and  vehicle  attitude  on  the  various  flight  param¬ 
eters.  This  was  done  through  the  use  of  a  computer  programmed  parametric  study. 

The  variation  in  retrovelocity  required  to  enter  at  a  given  flight  path  angle  is  shown  in 
Figure  2.  It  is  seen  to  increase  sharply  for  entry  angles  steeper  than  -0.5°.  Figure  3  shows 
the  optimum  angle  of  arrack,  and  Figure  4  shows  the  optimum  bank  angle  which  define  the 
optimum  entry  attitudes.  Under  these  conditions,  a  variety  of  (L/D)  ^  and  entry  angles 

were  run  on  the  computer  and  the  resulting  performance  is  shown  in  Figures  5  through  8. 
Notice  in  Figure  5  that  it  is  advantageous  to  force  the  velocity  vector  to  }  =  0°  before  it 
occurs  naturally,  though  for  increasing  L/D  one  may  wait  longer  and  take  advantage  of  the 
heading  change  being  achieved  without  using  thrust.  Note,  also,  that  the  higher  L/D  value 
causes  the  vehicle  to  penetrate  more  deeply  into  the  atmosphere  (Fig.  6),  yet  it  retains 
more  energy  (Fig.  7)  permitting  a  higher  plane  change  efficiency. 

In  considering  the  aerodynamic  maneuvering  portion  of  the  flight,  two  facts  become 
readily  apparent.  First  (Fig.  10),  steeper  optimum  bank  angles  result  from  low  lift- to- 
weight  ratios  and  high  L/D  values.  Second  (Fig.  12),  this  combination  results  in  greater 
plane  change  efficiency  by  reducing  drag  losses  and  directing  more  of  the  lift  force  toward 
plane  change  and  less  toward  keeping  the  vehicle  airborne.  Also  apparent  from  Figure  12 
is  the  fact  that  atmospheric  plane  change  efficiency  varies  nearly  linearly  with  increasing 
L/D. 


Since  the  higher  L/D  vehicles  penetrate  more  deeply  into  the  atmosphere,  it  might  be 
assumed  that  a  greater  velocity  impulse  would  be  required  to  exit  the  atmosphere.  Such  is 
not  the  case,  as  seen  in  Figure  13.  Despite  the  deeper  penetration,  substantially  less  impulse 
is  required  to  free  it  from  the  atmosphere  due  to  the  low  drag  coefficient  of  the  high  L/D 
vehicles.  Part  of  the  advantage  accumulated  by  the  higher  L/D  vehicles  is  erased  by  the 
velocity  requirement  for  orbital  reestablishment,  however,  as  seen  in  Figure  15.  While  a 
substantial  increase  does  exist,  it  is  by  no  means  excessive  and  does  not  greatly  affect  the 
overall  efficiency  which  appears  to  be  inherent  in  high  L/D  vehicles. 


The  advantage  gained  by  offsetting  the  thrust  is  seen  to  be  greatest  for  lower  L/D  values, 
diminishing  noticeably  as  L/D  increases  (Fig.  16).  This  advantage  is  taken  into  account  in 
the  succeeding  steps  which  define  the  optimum  flight  path  for  each  L/D  value. 


The  total  velocity  impulse  required  and  heading  change  achieved  for  all  phases  except  the 
atmospheric  maneuvering  (i.e.,  deorbit,  circularize  and  pullout,  exit  and  orbital  reestab¬ 
lishment),  including  offsetting  the  thrust  at  6  ,  were  compiled.  The  ratio  of  these  two 

quantities  (Ai/j/AV)  for  each  entry  condition  was  computed  and  plotted  in  Figure  19.  Note 
again  the  advantage  of  high  L/D  vehicles.  The  optimum  entry  angle  for  maximum  heading 
change  efficiency  is  obtained  from  this  figure.  The  locus  of  these  points  is  presented  in 
Figure  20.  This  optimum  maneuver  line  has  been  indicated  by  a  broken  line  on  all  preceding 
performance  curves. 
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Figure  19.  Heading  Change  Efficiency  for  Deorbit,  Pullout,  Exit  and  Reestablishment 
Including  Thrust  Offset  for  Various  L/D  Values 
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Figure  20.  Optimum  Entry  Angle  Versus  L/D 


Up  to  this  point  we  have  computed  heading  change  (A0)  and  not  orbital  inclination  change 
(Ai).  In  general,  change  of  heading  produces  part  inclination  change  and  part  orbital  node 
change,  the  ratio  depending  on  the  vehicle’s  location  with  respect  to  the  equator.  Exactly 
on  the  equator  heading  change  produces  only  inclination  change  while  90°  from  the  equator 
only  node  change  is  achieved.  From  the  consideration  of  the  interaction  of  the  various 
orbital  parameters,  it  becomes  inefficient  to  perform  heading  change  for  extended  time 
periods  as  shown  in  Figure  21.  For  example,  90  percent  efficiency  results  in  a  maximum 
turning  time  of  approximately  21  minutes.  Times  in  excess  of  this  result  in  substantial 
decreases  in  turning  efficiency,  as  seen  in  Figure  21.  Consequently,  it  is  necessary  to 
obtain  an  expression  for  the  time  spent  in  the  atmosphere,  and  then  the  total  time  con¬ 
sumed  by  a  given  plane  change.  The  time  for  the  turn  was  computed  from  an  expression 
for  the  thrusting,  banked  turn  derived  in  Appendix  III.  Since  the  reentry  and  exit  times 
can  be  obtained  directly  from  the  computer,  they  may  be  added  to  those  calculated  by 
Equation  (III- 16)  and  the  result  is  shown  in  Figure  22.  For  heading  changes  greater  than 
those  shown  in  Figure  22,  it  is  necessary  to  penetrate  more  deeply  into  the  atmosphere 
which  decreases  the  turning  time,  but  also  decreases  the  efficiency  (Figs.  6  and  19).  Since 
the  specific  purpose  of  this  report  is  to  define  the  most  efficient  maneuver,  the  above  cir¬ 
cumstances  have  not  been  considered  in  detail. 

One  now  reaches  the  point  where  it  is  possible  to  ascertain  whether  employment  of  aero¬ 
dynamic  maneuvering  and  its  attendant  complexity  is  really  worthwhile.  Figure  23  shows 
the  total  velocity  required  to  accomplish  a  given  plane  change  under  fully  optimized  condi¬ 
tions.  It  is  readily  apparent  that  the  savings  over  pure  propulsion  are  substantial.  In  the 
extreme,  2300  ft/sec  of  AV  allow  an  L/D  =  4  vehicle  to  attain  approximately  300  percent 
more  plane  change  than  for  pure  propulsion  in  orbit.  Without  question,  this  is  a  substantia! 
saving. 

Figure  24  shows  the  total  time  required  to  perform  a  given  plane  change  under  optimum 
conditions.  While  the  times  are  much  longer  than  for  pure  propulsion,  requiring  nearly  one 
full  orbit  for  17°  of  plane  change,  this  is  a  small  price  to  pay  for  the  propellant  savings 
attained. 
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Figure  21  Orbital  Inclination  Change  Efficiency  Versus  Time 


Figure  22.  Optimum  Turning  Time  in  the  Atmosphere  for  Various  L/D  Values 
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Figure  23.  Total  Velocity  Required  for  a  Given  Heading  Change  for  Various  L/D  Values 


Although  structural  limitations  have  not  been  considered,  it  is  interesting  to  note  what 
demands  are  made  on  the  structure  by  this  method.  Figure  25  depicts  the  peak,  equilib¬ 
rium,  stagnation  temperature  calculated  for  the  various  reentry  conditions  utilizing  the 
expression  (Ref.  5) 


\  =  [ -?!—  JT  f  1°  29  (37) 

\  1000  /  L  0.481  CD  \  1000/  -I  '  ' 

where  a  nose  radius  of  one  foot  and  an  emissivity  of  unity  have  been  assumed.  Note  that 
for  the  L/D  =  4  optimum  maneuver,  a  peak  temperature  of  4900°F  is  reached,  which  is 
not  beyond  the  forseeable  state  of  the  art  for  nose  cap  materials. 
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Figure  24.  Total  Time  to  Perform  a  Given  Heading  Change  for  Various  L/D  Values 
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IV.  CONCLUSIONS 


A  method  for  accomplishing  orbital  plane  change  utilizing  atmospheric  maneuvering 
with  propulsion  augmentation,  termed  the  “synergetic  or  aeropropulsive  method,”  has 
been  investigated.  The  performance  of  vehicles  having  L/D  values  from  1  to  4  was  com¬ 
pared  to  a  pure  propuision  plane  change. 

From  the  foregoing  analysis  the  following  conclusions  are  drawn: 

(1)  The  use  of  aerodynamic  maneuvering  in  performing  orbital  plane  change  results  in 
substantial  fuel  savings  over  plane  change  in  space  especially  when  one  follows  the  opti¬ 
mum  flight  path  and  flight  plan  described  herein  and  summarized  below. 

Deorbit  is  initiated  by  a  retrothrust  sufficient  in  magnitude  to  enter  the  atmosphere 

(by  a  Hohmann  ellipse)  at  the  optimum  flight  path  angle  shown  in  Figure  20.  The  vehicle 

is  flown  at  the  angle  of  attack  for  (L/D)  throughout  the  entire  synergetic  maneuver. 

max 

The  optimum  bank  angle  for  entry  and  exit  phase  is  shown  in  Figure  4.  When  the  vehicle 
reaches  the  optimum  altitude  (Fig.  6)  thrust  is  added  to  nearly  circularize  the  velocity 
and  the  bank  angle  increased  to  approximately  90°  (Fig.  11).  Thrust  is  then  added  to  main¬ 
tain  constant  speed  and  the  lift  vector  used  to  change  the  heading  in  the  “synergetic” 
turn.  After  the  desired  heading  change  has  been  accomplished,  a  velocity  impulse  (Fig.  15) 
is  added  to  return  to  the  original  orbit,  again  on  a  Hohmann  ellipse.  Finally,  velocity  is 
added  for  recircularization.  The  velocity  impulses  added  while  in  the  atmosphere  are  off¬ 
set  at  an  angle  from  the  original  vector  to  produce  higher  turning  efficiencies  (Fig.  16). 

(2)  Increasing  the  L/D  value  of  a  vehicle  substantially  improves  its  plane  change  effi¬ 
ciency,  thus  producing  greater  plane  change  for  each  pound  of  fuel  expended.  Specifically 
for  an  L/D  of  4  performance  over  pure  propulsion  in  orbit  was  improved  by  approximately 
a  factor  of  3. 

(3)  Although  maneuver  time  is  considerably  longer  than  for  pure  propulsion,  it  does  not 
become  a  constraining  factor  unless  plane  changes  greater  than  17°  are  desired.  In  this 
case,  it  would  be  necessary  to  perform  the  plane  change  in  stages,  with  perhaps  a  half 
orbit  being  required  between  stages,  if  the  efficiency  of  the  method  is  to  be  maintained. 

(4)  Although  the  temperatures  inherent  in  the  optimum  procedures  are  high  they  are 
within  the  foreseeable  state  of  the  art  for  advanced  materials. 

(5)  Unless  time  is  at  a  premium,  or  high  fuel  consumption  can  be  easily  tolerated,  the 
application  of  aerodynamic  maneuvering  to  orbital  plane  change  appears  worthy  of  con¬ 
sideration  for  future  space  missions. 
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APPENDIX  I 

DERIVATION  OF  AERODYNAMICS 


The  aerodynamics  used  in  this  report  are  based  on  the  work  of  Chapman  (Ref.  4).  Assum¬ 
ing  Newtonian  hypersonic  flow,  the  L/D  of  a  vehicle  may  be  written 


where 


.  2 

sin  a  cos  a 
b  +  sin3  a 


(1-1) 


'  b  = 


Cn  -  Cn 
umox  Dc 


(1-2)  ■ 


Cqo  =the  drog  coefficient  ot  a  =  0 


Cn  =  the  drog  coefficient  ot  a  =  90°. 
umax 

From  this  basic  definition,  several  quantities  may  be  obtained.  First,  it  is  desired  to  know 

the  variation  of  (L/D)  with  a  under  the  assumed  conditions.  At  a  d(L/D)/da  =  0 
v  '  'max  opt  opt  v  '  " 

and  hence,  from  Equation  (1-1) 

-(*) 


(I) 


=  2  cot  a  -  ton  a  — 


da 


3  sin  a  cos  a 
(b  +  sin3a  ) 


Setting  Equation  (1-3)  equal  to  zero,  an  equation  for  a  is  obtained, 


sina 


b  = 


opt 


2  cotza  .  -  I 


Then  (L/D)  is  defined  by 
v  '  'max  ' 


opt 


(t) 


mox 


[  2cotaont  “  ,onaopt  ]  ' 


(1-3) 


(1-4) 


(1-5) 


The  value  of  (WD)max  calculated  by  Equation  (1-5)  is  shown  in  Figure  I- 1 .  The  value  of  b 
corresponding  to  each  of  these  values  is  obtained  from  Equation  (1-4). 

To  obtain  the  C,  and  C^  values  for  each  of  the  (L/D)  values,  from  the  reference 

J  1_) 


max 

.  3 


Cn  -  Cr»  4-  (CD  CD  )  sin  a 

u  uo  '  umox  uo  • 

which,  using  Equation  (1-2)  becomes 


(1-6) 


Cq  =  (  -  )[b  +  sin3a|  , 

u  '  umox  uo  '  L  J 


(1-7) 
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Also,  from  Equation  (1-1),  combined  with  Equation  (1-7) 

Cj  =  (  Cn  -  Cn  )  sin2a  cosa 
*-  umax  uo 

From  the  assumption  of  Newtonian  flow, 


Hence 


and 


'max 


Cq  =  2  [  b  +  sin3 a  j 
CL  =  2  sin z a  cosa  . 


(1-8) 


(1-9) 


(I- 10) 

(1-11) 


Using  the  value  of  b  and  a  from  Equations  (1-4)  and  (1-5),  the  values  of  and 

at  (L/D)  are  obtained  from  Equations  (I- 10)  and  (I- 11).  These  values  are  tabulated  in 
max  ' 

Table  1-1  with  the  angle  of  attack  at  which  they  occur. 

Figure  1-2  has  been  placed  in  the  report  merely  for  the  convenience  of  the  reader.  It 
defines  the  drag  on  a  vehicle  having  the  L/D  value  shown  when  flying  at  the  indicated  alti¬ 
tude  and  at  local  circular  speed  over  the  equator. 


Figure  1-1.  (L/D)  Versus  Angle  of  Attack 
max 
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TABLE  1-1 


VEHICLE  CHARACTERISTICS  AT  (L/D) 

v  '  'max 


( o)max 

aopl 

(degrees) 

ci 

Lopt 

C°0Pt 

b 

0.33 

45 

0.707107 

2.12132 

0.707107 

0.5 

40.4 

0.639783 

1.279566 

0.367534 

1.0 

29.25 

0.41662 

0.41662 

0.09165 

1.5 

22.17 

0.26374 

0.17583 

0.03418 

2.0 

17.51 

0.17208 

0.08624 

0.01593 

2.5 

14.35 

0.11902 

0.04761 

0.00858 

3.0 

12.2 

0.0873 

0.0291 

0.0051125 

3.5 

10.6 

0.06652 

n  019006 

0.032785 

4.0 

9.35 

0.05209 

0.013022 

0.0022228 
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Figure  1-2.  Drag  Variation  with  Altitude  at  Circular  Velocity  for  Various  L/D  Values 
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APPENDIX  II 

PARAMETRIC  STUDY  OF  VEHICLE  REFERENCE  AREA  AND  MASS 


For  all  of  the  analyses,  a  constant  reference  area  and  initial  mass  have  been  assumed. 
This  assumption  is  justified  from  Figures  II- i  and  II- 2  which  show,  for  two  bank  angles, 
the  effect  of  mass  and  reference  area  variation  on  the  pullout  velocity.  This  parameter 
is  chosen  because  it  is  directly  related  to  the  energy  requirements  which  are  under  con¬ 
sideration.  It  is  seen  that  the  effect  of  the  variation  of  these  parameters  over  the  entire 
range  of  entry  angles  is  small  enough  to  be  neglected.  Thus  an  average  value  is  chosen 
and  held  constant  throughout  the  analysis. 
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Figure  II- 1.  Effect  of  Area  Variation  on  Pullout  Velocity 
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Figure  II-2,  Effect  of  Mass  Variation  on  Pullout  Velocity 
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APPENDIX  III 

TIME  REQUIRED  TO  PERFORM  A  PLANE  CHANGE 


Another  parameter  to  be  considered  in  this  evaluation  is  the  time  required  to  perform 
a  given  plane  change  utilizing  the  synergetic  method.  To  obtain  the  total  elapsed  time,  it 
is  necessary  to  evaluate  each  phase  independently  and  to  sum  the  results.  For  the  deorbit 
and  return  to  orbit  phases  the  equations  of  classical  celestial  mechanics  may  be  used  to 
obtain  the  ’'equired  time.  From  the  basic  equations  of  motion,  it  can  be  shown  that  the 
angle  swept  by  a  radius  vector  may  be  written 


ion  8  = 


V  siny 


7, 


0. 


(IH-1) 


where  all  values  are  instantaneous  and  0  is  measured  from  the  orbit  perigee.  By  using 
basic  definitions,  Equation  (111- 1)  may  be  written 


tan  6 


w2 

V  siny  cosy 

T72  2  ,, 2  ' 

V  cos  y  -  Vc 


(HI-2) 


If  Equation  (III- 2)  is  chosen  to  represent  the  deorbit  ellipse,  following  the  firing  of  the 
retrograde  impulse  while  in  a  given  circular  orbit,  then  the  angle  swept  from  retrofire 
to  reentry  may  be  rewritten,  in  terms  of  the  reentry  and  original  orbit  parameters  as 


ton  6 


l  Ven  \2 

l^T)  Sm>'en  cos7en 


(111-3) 


'i 


en 


»  ven 

hr  005 


7  )  -- 
fen  J  r 


f±_ 

en 


The  eccentricity  of  the  reentry  orbit  can  be  obtained  from  the  following  equation  ex¬ 
pressed  in  terms  of  the  radii  at  the  apogee  and  perigee; 


2  r. 


rA  +  rp 

Using  the  definition  of  the  semi- major  axis  of  an  ellipse, 


(111-4) 


2o  =  rp  +  rA 

the  relation 


and  the  definition  of  circular  velocity 


(111-5) 


(111- 6) 


v.  =  4- 


(111-7) 
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Equation  (III— 4)  may  be  written,  afteze-  considerable  manipulation  as 


(II1-8) 


To  obtain  the  time  requited  for  the  vehicle  to  start  reentry  after  the  retrograde  impulse, 
Kepler’s  law  of  equal  areasswept  by  a  radius  vector  in  equal  times  yields 


or 


d_A 

dl 


=  Constant 


(HI-9) 


dl 


d  B 


where  0  is  the  angle  swept, 


Since,  in  general,  the  radius  r  is  given  by 


I  +  €  COS  8 


(III- 10) 


(111-11) 


Equation  (111-10)  may  be  integrated,  ior  an  ellipse  (0  <  e  <  1)  to  yield 


to  “  1| 


g^3  f  sinf 

-I  )  ( lie  cosQ  } 


( I  —  c  ) 

/T7 


By  virtue  of  the  fact  that  on  a  circular  orbit  the  values  of  a  parameter  at  all  points  are  the 
same,  the  point  of  retrofit  will  be  taken  as  the  apogee  which,  in  fact,  it  is  for  the  trans¬ 
fer  ellipse.  Thus  0:  =  tt  bj definition  and  Equation  (III- 12)  becomes,  since  q  =  0,  after 
some  manipulation 


en 


ren  ven  C0S7e> 

Vc2  (1+  €.  ) 
CA  2 


vT 


tan 


\  (|  +  e# 


where 


(111-13) 


0gn  is  defined  by  Equation (m_ 3)  and  e2  by  Equation  (111-8). 

By  treating  the  exit  maieuver  as  the  reverse  of  reentry,  Equations  (111-3),  (1II-8)  and 
(111-13)  may  be  used  to  calculate  ttae  exit  conditions  and  the  apogee  values  of  the  param¬ 
eters  for  the  completion  j[ the  maneuver,  The  time  required  to  enter  the  atmosphere  and 
pullout  and  exit  the  atmosphere  are  obtained  directly  from  the  computer  runs. 
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The  time  required  to  perform  the  atmospheric  maneuver  may  be  obtained  from  Equations 
(26),  (24),  (15),  (19)  and  (23).  Integration  of  Equation  (23)  determines  the  mass  ratio 


J  n 


VAi// 


Vex  si 


+ 1 


Assuming  m  is  constant,  then 


At 


Am  m2  -  m. 


man 


(III- 14) 


(III- 15) 


Thus,  utilizing  Equations  (III- 15),  (III- 14),  (15)  and  (19),  the  maneuver  time  may  be  written 
as 


At 


m,  Vex  cosA 


man 


[  I  -  exp  (  - 


VA  \(/ 


vex  s 


■)  i.  (in-IQ 


+i 


The  times  calculated  by  Equation  (III- 1 6),  added  to  the  time  required  to  enter  and  exit  the 
atmosphere  are  shown  in  Figure  13.  A  useful  engineering  approximation  for  small  angles 
is 


(2.35) 


minutes 


_L_ 

W 


(III- 1 7) 
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